Resonance width distribution for high-dimensional random media 
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We study the distribution of resonance widths V(T) for three-dimensional (3D) random scattering 
media and analyze how it changes as a function of the randomness strength. We are able to identify 
in V(T) the system-inherent fingerprints of the metallic, localized, and critical regimes. Based on the 
properties of resonance widths, we also suggest a new criterion for determining and analyzing the 
metal-insulator transition. Our theoretical predictions are verified numerically for the prototypical 
3D tight-binding Anderson model. 
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I. INTRODUCTION 

Quantum mechanical scattering in systems with com- 
plex internal dynamics has been a subject of intensive 
research activity for a number of years. The inter- 
est was motivated by various areas of physics, rang- 
ing from nuclear^, atomioS. and molecular— physics, 
to mesoscopics 4 , quantum chaos 5 * 6 ", and classical wave 
scattering 7 . Recently, the interest in this subject was 
renewed due to technological developments in quantum 
optics associated with the construction of new type of 
lasers^ and the experimental investigation of atoms in 
optical lattices^. 

The most fundamental object which characterizes the 
process of quantum scattering is the scattering matrix 
S, where S relates the amplitudes of waves that en- 
ter and leave a scattering region. Of great interest 
are the statistical properties of the poles of the S ma- 
trix. They determine the conductance fluctuations of a 
quantum dot in the Coulomb blockade regimeii or the 
current relaxation 1 ^. The poles of the S matrix are 
related to resonance states occurring at complex ener- 
gies £ n — E n — ^r„, where E n is the position and T n 
the width of the resonance. Resonances correspond to 
"eigenstates" of the open system that decay in time due 
to the coupling to the "outside world" . 

For chaotic systems Random Matrix Theory (RMT) 
is applicable and the distributions of resonance widths 
V(T) is known. A review can be found in Ref. (see 
also Ref. fl4j) . As the disorder increases, the system be- 
comes diffusive and the deviations from RMT increase 
drastically. For low dimensional random systems in the 
metallic regime the distribution of resonances V(T) was 
found recently 15-1617 . For the strongly disordered limit, 
where localization dominates, V(Y) was investigated by 
various groupsi&isLSii as well. At the same time an at- 
tempt to understand systems at critical conditions was 
undertaken in Refs. |^ and |2j The latter deals with 
random systems of higher dimensions, the most promi- 
nent of which is the three-dimensional (3D) Anderson 



model. It undergoes a Metal-Insulator Transition (MIT) 
with increasing strength of disorde»2£. 

In this paper we extend our previous analysis on the 
3D Anderson model 2 ^ and study the distribution of reso- 
nance widths as we change the disorder strength. Based 
on the analysis of V(T) we propose a new method to 
locate the MIT. The paper is organized as follows: In 
Section [H] the 3D Anderson model and the scattering 
formalism are introduced. In Sec. IHII we discuss the con- 
sequences of localization in the distribution of resonance 
widths in the diffusive and localized regimes as well as at 
the MIT and show the numerical results supporting our 
arguments. In Sec. IIVI we investigate a new method for 
determining and analyzing the emergence of the MIT and 
propose a scaling theory near the critical point. Finally, 
our conclusions are given in Sec. 



II. THE 3D ANDERSON MODEL AND THE 
SCATTERING SETUP 

The Anderson model with diagonal disorder on a 3D 
cubic lattice is described by the tight-binding Hamilto- 
nian (TBH) 



H = J2\n)W n (n\+ £ l n X m l 



(1) 



(n,: 



where n = (n x ,n y ,n z ) labels all the N — i 3 sites of 
the cubic lattice, while the second sum is taken over all 
nearest-neighbor pairs (n, m) on the lattice. The on-site 
potential W n for 1 < n x , n y , n z < L is independently and 
identically distributed with probability V(W n ). We use 
three different distributions for the random potential: (a) 
a box distribution, i.e., the W n are uniformly distributed 
on the interval [-W/2, W/2]; (b) a Gaussian distribution 
with zero mean and variance W 2 /12; and (c) a Cauchy 
distribution V{W n ) = W/-k(W 2 + W 2 ). For the system 
defined by Eq. JTJ the MIT for E ~ occurs for W = W c 
with (a) W c ~ 16.5, (b) W c ~ 21.3, and (c) W c ~ 4.26 
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FIG. 1: Scattering setup. The sample is a cubic lattice of 
linear length L. To each of the M = L 2 sites of the layer 
Tlx = 1 semi-infinite single mode leads are attached. 



(see Ref.H3). Then, for W < W c {W > W c ) the system 
is in the metallic (insulating) regime. 

We turn the isolated system to a scattering one by 
attaching M = L 2 semi-infinite single mode leads to each 
site of the layer n x = 1, as depicted in Fig. ^ Each lead 
is described by a one-dimensional semi-infinite TBH 

— oo 

H M = ^2(\n><n + l\ + \n+l><n\) . (2) 

n=l 

Using standard methods^ one can write the scattering 
matrix in the form 21,22 

S(E) = l-2ism(k)W T (E-H eS )~ 1 W , (3) 

where 1 is the M x M unit matrix, k — arccos(i^/2) is 
the wave vector supported in the leads, and H c s is an 
effective non-hermitian Hamiltonian given by 

n cS = H -e lk WW T . (4) 

Here, W is a N x M matrix that specifies at which site of 
the sample we attach the leads. Its elements are equal to 
zero or y/w, with < \/w < 1, where w is the coupling 
strength. Below, unless stated otherwise, we will always 
consider the case w = 1. Moreover, since arccos(£'/2) 
changes only slightly in the center of the band, we set 
E = and neglect the energy dependence of Ti c s- The 
poles of the S matrix are then equal to the complex zeros 
of 

det[£ - H eS ] = 0. (5) 

^From Eqs. 10 and J^J it is clear that the formation 
of resonances is closely related to the dynamics in the 
scattering region, governed by Hq. 

In order to investigate the distributions of resonance 
widths we used samples with L = 20 as a maximum size. 
For better statistics a considerable number of different 
disorder realizations was considered. In all cases we had 
at least 10 000 data for statistical processing. 



III. DISTRIBUTION OF RESONANCE 
WIDTHS: RESULTS AND DISCUSSION 

A. Metallic Regime 

When the disorder strength W is smaller than W c , but 
still large enough so that the mean free path is smaller 
than the system size, the system is in the metallic regime. 

Recently, a lot of research activity was devoted to 
the understanding of the statistical properties of various 
physical quantities (such as conductance, local density of 
states, current relaxation times) in finite-size random sys- 
tems in the metallic regime. The outcome of these stud- 
ies indicated that the tails of these distribution functions 
show large deviations from the universal Random Matrix 
Theory (RMT) results, expected to be valid 2 ^ in the limit 
of infinite dimensionless conductance g = Txh/A = DL. 
Here, Txh ~ D/L 2 is the typical inverse time (Thouless 
time) that an excitation needs to diffuse (with diffusion 
coefficient D) in order to reach the boundary of a sys- 
tem, with linear size L, and A ~ 1/L 3 is the mean level 
spacing. 

The origin of these deviations was found to be related 
to the existence of eigenstates which are unusually lo- 
calized around a center of localization. These states are 
precursors of the Anderson localization and were termed 
prelocalized states^SiSL 2 ^ 2 ^. In 3D conductors they have 
sharp amplitude peaks on the top of a homogeneous 
background^Si. 

We start our analysis by investigating the effects of pre- 
localized states in the distribution of resonance widths. 
It is natural to expect that these states with localization 
centers at the bulk of the sample are affected only weakly 
when opening the system at the boundaries. Therefore, 
prelocalized states decay very slowly to the continuum 
leading us to the conclusion that the corresponding res- 
onance widths (inverse lifetime) T are smaller than the 
mean level spacing A. Hence, assuming the validity of 
standard first order perturbation theory (that can be ap- 
plied if the coupling of the sample to the leads is weak, 
w <C 1) we get 

^ = <¥|WtW|¥) oc Yl \*(n)\ 2 ~ L 2 \*(L)\ 2 , 

nGboundary 

where ^(L)! 2 is the wavefunction intensity of a pre- 
localized state at the boundary of the sample. At the 
same time, the distribution of wavefunction components 
at the boundary was found to be2£ 

P(0)~exp[-Ci]n 3 (0)] , (7) 

with 0- 1 = L^f(L) and C x oc g. Using Eq. J7J) together 
with Eq. © we obtain 

7>(l/r)~exp[-C 2 ln 3 (l/r)] , (8) 

where C2 oc g. 
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ln 3 (l/T) 

FIG. 2: The distribution of resonance widths, plotted as 
P(l/r) vs ln 3 (l/r), for T < A in the diffusive regime. L = 16 
and W = 10, 12, and 14 (from left to right). The dashed lines 
are (shifted) linear fittings to the distributions. 




FIG. 3: The resonance width distribution V(T) for L = 16 
and various scattering configurations (from up to down): M = 
L 2 leads attached to the boundary, one lead attached to the 
boundary, and one lead attached to a site in the bulk of the 
sample. The dashed lines are the corresponding theoretical 
predictions for F > Fxh given by Eqs. 11UI12II . see main text. 



As can be seen from Fig. [2] the prediction of Eq. JSJ, 
obtained using perturbation theory, holds even for strong 
coupling. Indeed, the reported data for the 3D Ander- 
son model in the diffusive regime, plotted as \nV(l/r) 
vs ln 3 (l/r), shows a linear behaviour. This comes as a 
surprise since in Fig. |2 we have considered perfect cou- 
pling, w = 1. Differences in slope correspond to different 
dimensionless conductances g induced by the varying val- 
ues of disorder strengths. 

Next, we turn to the analysis of V(T) for T > Txh 
A. In order to go on we need to recall that the inverse of T 
represents the quantum lifetime of a particle in a resonant 
state escaping into the leads. Moreover, we assume that 
the particles are uniformly distributed inside the sample 
and spread until they reach the boundary where they 
are absorbed. Then, we can associate the corresponding 
lifetimes with the time t R ~ l/r# a particle needs to 
reach the boundaries, when starting a distance R away. 
The relative number of states that require a time t < 
t R in order to reach the boundaries (or equivalently the 
number of states with T > T R ) is 

Vint(T R ) = J™ V(Y)dT ~ , ( 9) 

where V(t R ) ~ L 3 — (L — R) 3 is the volume populated 
by all particles with lifetimes t < t R . 
Assuming now diffusive spreading, 

R 2 = D ■ t R , 

we get from Eq. © (to leading order with respect to 

r Th /r) 

— v 3/2 2 y*2 

T -f) = (?) (10) 

where we refer to the rescaled variable T = T/A. Equa- 
tion l|lUfl is valid as long as the leads are attached to the 



boundary of the sample. 

Here, it is interesting to point out that a different way 
of opening the system might lead to a different power 
law behavior for V(T). Such a situation can be realized 
if instead of opening the system at the boundaries we 
attach one lead somewhere in the sample. In such a case 
we have 

leading to 

p(r)~r~ 5 / 2 , (ii) 

where we used V(t R ) ~ R 3 . The above results are correct 
for any number of leads M such that the ratio M/L 3 
scales as 1/L 3 . 

If, on the other hand, we attach a single lead to the 
boundary of the sample we obtain 

-p(r) - T~ 2 . (12) 

This behavior is due to the fact that the decay takes place 
at the surface, leading to a situation similar to that of a 
2D systemii. 

In Fig. |3]we present numerical data for the 3D Ander- 
son model in the metallic regime. We verify the validity 
of the theoretical predictions given by Eqs. l|lUll2fl by 
using various configurations: M — L 2 leads attached to 
the boundary, one lead attached to the boundary, and 
one lead attached to a site in the bulk of the sample, 
respectively. We observe a good agreement with the ex- 
pected behavior in all cases. 

B. Localized Regime 

When the disorder strength W is larger than W c the 
system is in the localized regime. In this regime the eigen- 
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functions are exponentially localized in space and, as a 
consequence, transmission is inhibited and the system 
behaves as an insulator. 

Various groupsi&iSiSii had investigated the resonance 
width distribution of low dimensional random media in 
the localized regime during the last years. In the region of 
exponentially narrow resonances r < Tq = e~xjp(—2L/l OQ ) 
the distribution was found to be log-normal, i.e., 



V{T) ~ exp 



in 2 (r) 



r<r n 



(13) 



This result is analogous to the conductance distribution 
of localized systems. Equation essentially relies on 
two assumptions: first, that eigenfunction components 
are randomly distributed with no long-range correlations; 
and second, that they are exponentially localized with a 
normal distribution of localization lengths. 

It is reasonable to assume that the same arguments 
leading to Eq. (|13[1 applies as well for high-dimensional 
random media like the 3D Anderson model in the local- 
ized regime. Indeed, our numerical results reported in 
Fig. Ufa) show good agreement with the theoretical ex- 
pectation l|13|) which supports our assumption. 

On the other hand, we found [see Fig. E£b)] that the 
long tails of the distribution behave as 



V(T) 



T < T < 1/L 



(14) 



Equation (|14fl can be easily understood when employ- 
ing Eq. lO- The new ingredient is that wavefunctions 

are exponentially localized: |^(f)| ~ & 3 ^ 2 exp(— r//oo). 
Using simple perturbation arguments^ [see Eq. ©], i.e. 
T ~ |*(r)| 2 , we obtain 

i? 3 ~4ln 3 (4r) ■ 

By inserting this into Eq. 10 we get Eq. iQ , to leading 
order with respect to loo/L. 

The region of large T values is essentially determined 
by the coupling to the continuum, so it should be model- 
dependent. Nevertheless, it is reasonable to assume that 
the number of resonances involved is constant, of order 
l OQl and therefore the extreme tail of the distribution 
should subside at large L at the rate ~ loo/L. 

Let us finally note that in the thermodynamic limit 
L — > oo the probability of finding an eigenstate at any 
finite distance from the boundary is equal to zero. Thus 
the distribution of the resonance widths in this case ap- 
proaches a delta function centered at zero. 



C. Criticality 

The MIT, where W = W c , is characterized by several 
critical properties: the level statistics acquires a scale- 
independent formSSiSSiSiiS while the eigenfunctions show 
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FIG. 4: V(T) in the localized regime for various combina- 
tions of W and L in the range T < 1. The log- normal de- 
cay is highlighted by Gaussian fits (full curves) whose maxi- 
mum decreases with increasing strength of disorder and also 
shifts towards smaller values of V. When keeping the ratio 
loo/L pa 0.136 fixed, similar distributions were obtained for 
different combinations of L and W (filled circles and open 
squares), (b) For T > 1 the power-law decay V(T) ~ 1/T 
is observed (dashed line) which becomes more robust for in- 
creasing strength of disorder. 



strong fluctuations on all length scales and obey multi- 
fractal distributiona 26 i 27 i 33 i 34 i 35 . 

In Ref. |52 it was found that V(T) follows a new uni- 
versal distribution, i.e., independent of the microscopic 
details of the random potential and number of attached 
leads. Specifically, it decays asymptotically with the 
power 



V(f) ~ 5 l/3f-(l+l/3) 



(15) 



which is different from those found for chaotic, metallic, 
or localized systems (see Ref. U| Eq. (fTUI , and Eq. 1(11)1 , 
respectively) . 

One can relate the power-law decay (|15f) to the anoma- 
lous diffusion at the MIT. Indeed, at the MIT the con- 
ductance of a 3D disordered sample has a finite value 
<7 C ~ 1. Approaching the MIT from the metallic regime 
one has g ~ £r/A, where Et = D/R 2 is the Thouless 
energy, D is the diffusion coefficient, and A ~ 1/i? 3 is 
the mean level spacing of a sample with linear size R. 
This yields D ~ g c /R at W c . Taking into account that 
D = R 2 /tn, we get for the spreading of an excitation at 
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FIG. 5: Universal behavior of V(T) at the MIT [reported 
here as 7 3 (ln(r))] for various sample sizes L and potential 
distributions. 



the MIT 



9c ■ t R 



Then, straightforward application of Eq. © leads to 
Eq. (BJ. 

In Figs. |5] and we report some numerical results for 
the 3D Anderson model at the MIT. Figure shows 
the distribution of the logarithm of the rescaled reso- 
nance widths , P(ln(r)) for the three different distribu- 
tions V(W n ) of the random potential and for various 
sample sizes L. The body of the distribution function 
in all cases coincides and does not change its shape or 
width. Of course, the far tail of this universal distribu- 
tion develops better with increasing L. The sharp peak 
appearing at the right is an artifact of our choice to ne- 
glect the energy dependence of H c s- We thus confirm 
that at the MIT the distribution of rescaled resonances 
is indeed scale-invariant independent of the microscopic 
details of the potential. 

/From Fig. EJa) an inverse power law 'Pint(r) ~ T~ a 
is evident. The best fit to the numerical data yields 
a = 0.333 ± 0.005 in accordance with Eq. I|15l) . Here, 
the case of perfect coupling (w = 1) has been consid- 
ered. Different coupling strengths are going to affect this 
behavior as can be seen from Fig. ffib). 



IV. A SCALING THEORY FOR THE 
RESONANCES WIDTHS 

In the original proposal of the scaling theory 
of localization, the conductance g is the relevant 
parameter 23,36 . A manifestation of this statement is seen 
in Eqs. (|8I10I13I14I15|) where V° nt = Vi n t(^o) is propor- 
tional to the conductance g. It is therefore natural to 
expect that V® ut will follow a scaling behavior for finite L 
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FIG. 6: (a) The integrated distribution 7-mt(r) in the case of 
perfect coupling, w = 1, for various sample sizes L and poten- 
tial distributions. The dashed line is the theoretical prediction 
Ptot(f) ~ f" ' 333 , see Eq. HHJ. (b) P int (f) for a box distri- 
bution, L = 10, and different coupling strengths: w = 0.001, 
0.01, and 0.5 from left to right. 



(and for some To ~ 1) that is similar to the one obeyed 
by the conductance g. The following scaling hypothesis 
was therefore postulated in Ref. 22: 



V? nt (W,L) = f(L/l 00 (W)) 



(16) 



In the insulating phase (W > W c ) the conductance of 
a sample with length L behaves as g(L) ~ exp(— L/Zoo) 
due to the exponential localization of the eigenstates, and 
therefore we have g{L\) < g(L 2 ) for L\ > L 2 . Based 
on Eq. I|15(l we expect the same behavior for V^ nt ; i.e., 
for every finite L\ > L 2 we must have V^ nt (W, L\) < 
Vi nt (W, L-z)- On the other hand, in the metallic regime 
(W < W c ) we have that g(L) = DL and therefore we ex- 
pect from Eq. dJ) Vf nt {W,L x ) > V? nt {W,L 2 ). Thus, the 
critical point is the one at which the size effect changes 
its sign, or in other words, the point where all curves 
V-* nt (W,L) for various L cross. One can reformulate the 
last statement by saying that in the thermodynamic limit 
L — > oo at W — W c the number of resonances with width 
larger than the mean level spacing goes to a constant. 

In Fig. El we show the evolution of Vi nt (W) for dif- 
ferent values of L using the box distribution. From 
this analysis the critical disorder strength W — W c — 
16.5 ± 0.5 was determined in agreement with other 
calculations 24 . A further verification of the scaling hy- 
pothesis l|16fl is shown in Fig. [H] where the same data 
are reported as a function of the scaling ratio L/l^. 
Note that all points collapse on two separate branches 
for W < W c and W > W c . 
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FIG. 7: Pi nt (W,L) as a function of W for different system 
sizes L provides a means to determine the critical point W c 
of the MIT (vertical line at W = 16.5). 
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confirmed for various system sizes L and disorder strengths 
W using the box distribution. 



V. CONCLUSIONS 



We have studied the properties of the resonance width 
distribution in various regimes for the 3D Anderson 
model. 

In the metallic regime we obtained the forms of V(T) 
for small and large T and show that they are determined 
by the underlying diffusive dynamics and by the exis- 
tence of prelocalized states. For the localized regime we 
also explored the limits of small and large T. In the first 
limit we found that V(T) shows a log-normal behavior 
while in the latter the distribution is power-law like. At 
the MIT we show that V(T), with f = T/A, has a uni- 
versal form, i.e., independent of the microscopic details 
of the random potential and number of attached leads. 
Specifically, it decays asymptotically with a power which 
is different from those found in the diffusive and local- 
ized regimes. In addition, based on resonance widths, we 
suggested a new method for determining and analyzing 
the emergence of the MIT and propose a scaling theory 
near the critical point. 
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